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ASYMPTOTIC ANALYSIS OF LOSS PROBABILITIES 
IN GI/M/m/n QUEUEING SYSTEMS AS n INCREASES 

TO INFINITY 

VYACHESLAV M. ABRAMOV 

Abstract. The paper studies asymptotic behavior of the loss 
probabihty for the GI/M/m/n queueing system as n increases 
to infinity. The approach of the paper is based on appUcations of 
classic results of Takacs (1967) and the Tauberian theorem with 
remainder of Postnikov (1979-1980) associated with the recurrence 
relation of convolution type. The main result of the paper is associ- 
ated with asymptotic behavior of the loss probability. Specifically 
it is shown that in some cases (precisely described in the paper) 
where the load of the system approaches 1 from the left and n 
increases to infinity, the loss probability of the GI/M/m/n queue 
becomes asymptotically independent of the parameter m. 
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1. Introduction 

It is well-known that queueing systems with many servers are well mod- 
els for communication systems. Study of queues with many servers and 
especially analysis of the loss probability have a long history going back 
to the works of Erlang (see [7J), who in 1917 first gave fundamental re- 
sults for Markovian queueing systems, and to the works of Palm, Pol- 
laczek and other researchers (e.g. [I5], [16], [H], [I9], [22]), who then 
developed the Erlang's results to non-Markovian systems. Since then 
these results have been developed in a large number of investigations, 
motivated by growing development of modern telecommunication sys- 
tems. Nowadays the theory of loss queueing theory is very rich. There 
is a number of different directions of the theory including management 
and control, redundancy, analysis of retrials, impatient customers and 
so on. 

In the present paper we study the Gl jMjmjn queueing system in 
which parameter n, the number of possible waiting places, is a large 
value. This assumption is typical for real telecommunication systems. 
Under this assumption the paper studies asymptotic behavior of the 
loss probability, where the most significant result seems to be related 
to the case where the load parameter approaches 1 from the left. The 
case of a heavy load parameter is the most interesting in practice. In 
the pre-project study and system design stage an engineer is especially 
interested to know about the behavior of the loss probability in heavy 
loaded systems. 

We consider the GI / M /m/n queue, where m > 1 is the number of 
servers, and n > is the admissible queue-length. Let A{x) denote 
the probability distribution function of an interarrival time, and let 
A be the reciprocal of the expected interarrival time. For real s > 
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we denote a{s) = e~'^^dA{x). The parameter of the service time 
distribution is denoted /i, and the load of the system is g = X/ {mjj). (In 
Theorem 13.21 the parameter q is assumed to depend on n. However we 
do not write this dependence exphcitly, assuming that this dependence 
is clear from the formulation of the aforementioned theorem.) 

In the case of GI /M/l/n queueing system for the stationary loss 
probability p„ we have the representation (Abramov ^) 

Pn = — , (1.1) 

where the generating function n(z) of tTj, j = 0, 1, ... is the following: 
n(.) = Vvr,.^ = -^^/i^i^, |.|<a, (1.2) 

a is the least in absolute value root of the functional equation z = 
a{fi — fiz) (the variable z is assumed to be real). It is well-known (e.g. 
Takacs ^23j), that cr belongs to the open interval (0,1) if g < 1, and it 
is equal to 1 otherwise. Note, that another representation for the loss 
probability pn is given in Miyazawa [Hj . 

The value 7r„ has the following meaning. This is the expected number 
of arrivals up to the first loss of a customer arriving to the stationary 
system. 7r„ satisfies the recurrence relation of convolution type 

n 

TTn = y^JiTln-i+l, n = 0, 1, . . . , (1.3) 
1=0 

where the initial value is ttq = 1. Specifically, 

= r e-^'^^^^dAix). (1.4) 
Jo 

As n —>■ oo, the asymptotic behavior of the general recurrence relation 
of convolution type, 

n 

Qn = fiQn-i+1, (l-^) 

1=0 
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where /o > 0, fi > (i > 1) and /o + /i + ••• = 1, has been origi- 
nally studied by Takacs [2^, p. 22 and then developed by Postnikov 
[17j . Section 25. For the readers' convenience, some of these results, 
necessary for the purpose of the paper, are collected in the Appendix. 

By exploiting (11. 3p . Abramov [2J studied the asymptotic behavior of 
the loss probability p„, as ^ oo. The analysis of Abramov [2] is based 
on the mentioned results of Takacs [SJ and Postnikov . For other 
applications of the mentioned results of Takacs [23] and Postnikov [17] 
see also Abramov [1], [3], [3]. 

The asymptotic analysis of GI / M/m/n queueing system, as n ^ 
oo, is a much more difficult problem than the same problem for the 
Gl jMjXjn queue. Recently, Choi et al [9] and Kim and Choi [12] 
obtained some new results related to the GI / M/m/n and Gl"^ /M/m/n 
queues. In Choi et al [9j the exact estimation for the convergence rate of 
the stationary GI / M /m/n queue-length distribution to the stationary 
queue-length distribution of the GI/M/m queueing system, as n ^ 
cxD, is obtained. In Kim and Choi [12] detailed analysis of the loss 
probability of the GI-^ /M/m/n is provided. The analysis of these two 
aforementioned papers is based on the development of the earlier results 
of Choi and Kim [8J in a nontrivial fashion. 

The analysis of Choi and Kim |8j and Choi et al [IQJ in turn uses the 
deepen theory of analytic functions, including an untraditional result 
of the theory of power series given by a theorem of Wiener. 

Ramalhoto and Gomez-Corral [18] discuss retrials in M/M/r/d loss 
queues and present an appropriate decomposition formulae for losses 
and delay in those queueing systems. In the case of Markovian system 
with many servers the results of Ramalhoto and Gomez-Corral [18] 
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are useful for analysis of the effect of retrials in asymptotic analysis of 
losses. 

The present paper provides the asymptotic analysis of the loss proba- 
bility of the GI / M /m/n queue as ^ oo, by reduction of the sequence 
'^m^,Q,n to above representation (11.51) . and then estimates the loss proba- 
bility Prn,n = '^/''^^mo n usiug the rcsults of Takacs [21] and Postnikov 
p!7] given in the Appendix. (The precise definition of the sequence 
^mi,n is given later.) 

This is the same idea as in the earlier paper of Abramov [2] re- 
lated to the GI / M /1/n queue, however it is necessary to underline the 
following. Whereas in the case of GI / M/l/n the reduction to ( 11. 3p is 
straightforward, and the representation of the probabilities r„ and their 
generating function is very simple, the reduction to the recurrence re- 
lation of (II. 5p in the case of GI / M /m/n queue, m > 1, is not obvious, 
and representation for probabilities rk^m_k^n and tq^^j and associated 
generating function is more difficult. (The aforementioned probabilities 
fk,m-k,n and ro,m,j are defined later.) Furthermore, the sequence T^mfij^ 
j < n, is structured as schema with the series classes, and, as j = n, the 
value pm,n = ^/'^mfl,n coiucides with the desired loss probability. For 
this reason the class of asymptotic results, that we could obtain here 
for GI /M/m/n queue, is poorer than that for the GI / M /1/n queue in 
Abramov [2]. 

Our approach has the following two essential advantages compared 
to the pure analytical approaches of Choi et al [lO], Kim and Choi [T2] . 
Choi et al [9j, Simonot [20] : 

• The problem reduces to the known classic results (Theorem of 
Takacs [24] and Tauberian Theorem of Postnikov pJZj), permitting us 
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to substantially diminish the cumbersome algebraic calculations and 
clearly understand the results. 

• Along with standard asymptotic results having a quantitative fea- 
ture we also prove one interesting property related to the case of n 
increasing to infinity and the load approaching 1 from the left. (For 
the more precise assumptions see formulation of Theorem 13. 2[ ) Specif- 
ically it is proved that the obtained asymptotic representation is the 
same for all m > 1, i.e. it coincides with asymptotic representation 
obtained earlier for the GI / M/l/n queueing system in Abramov [2J. 
As n increases to infinity, this asymptotic property remains true for all 
fixed ^ > 1. 

The conditions of Theorem 13.21 and the first two cases of Theorem 
13.11 all fall into the domain of heavy traffic theory (e.g. Borovkov [6J, 
Whitt [25], [26]). Although the aforementioned recent results of Whitt 
|25j . [26] are related to more general models, they however do not cover 
the results of this paper. 

Our approach is based on asymptotic analysis of relations (13. 4p and 
(13.51) which is based on asymptotic representation for the root of equa- 
tion z = a{fim—fimz) (see formulation of Theorem 13. ip as g approaches 
1, presented in the book of Subhankulov [21]. (Chapter 9 of this book 
is devoted to application of Tauberian theorems to a specific moving 
server problem arising in Operations Research.) 

The rest of the paper is organized as follows. In Section 2 we give 
some heuristic arguments preparing the reader to the results of the 
paper. There are two theorems presenting the main results in Section 
3. In Section 4 we derive the recurrence relation of convolution type for 
the loss probability. These recurrence relations are then used to prove 
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Theorem 13.11 of the paper. Theorem 13.11 is proved in Section 5 and 
Section 6. In Section 7 we study the behavior of the loss probabihty 
as the load approaches 1 from the left. In Section 8 numerical example 
supporting the theory is provided. The Appendix contains auxiliary 
results necessary for the purpose of the paper: the theorem of Takacs 
|24j and the Tauberian theorem of Postnikov |17j . 

2. Some heuristic arguments 

The GI /M/m/n queueing system, m > 1, is more complicated than its 
analog with one server. Letting n be infinity, discuss first the GJ/M/l 
and GI /M/m queueing systems. 

It is well-known that the stationary queue-length distribution of 
GI/M/1 queue (immediately before arrival of a customer with large 
order number) is geometric. The same (customer-stationary) queue- 
length distribution of GI /M/m queue, provided that immediately be- 
fore arrival at least m — 1 servers are occupied, is geometric as well. 
Thus, a typical behavior of the queue- length processes of GI/M/1 and 
GI /M/m queues is similar, if we assume additionally that a customer 
arriving into GI /M/m queueing system finds at least m — 1 servers 
busy (see Kleinrock [ISj). 

The similar situation holds in the case of the GI / M/ 1/n and GI / M/ m/n 
queues. Specifically, in the case of the GI / M /1/n queue the station- 
ary loss probability satisfies f ll.ip - fll.3l) . In the case of the GI /M/m/n 
queue the conditional stationary loss probability provided that upon ar- 
rival at least m — 1 servers are busy satisfies (11. 3p as well. In the case 
of the GI /M/m/n queue the only difference is that, the value fi in 
(11.21) and (11.41) should be replaced with fim, and a should be the least 
in absolute value root of the equation z = a{fj,m — fj,mz) rather than 
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of the equation z = a{fi — jiz). In the sequel, the least root of the 
functional equation z = a{jj,m — fimz) is denoted am- 

Let us now discuss the stationary probabilities of the GI /M/m/n 
system, m > 1, without the condition above. It is clear that eliminating 
the condition above proportionally changes the stationary probabilities 
P{ arriving customer meets m+j — 1 customers in the system }, j > 0. 
That is, the loss probability is changed proportionally as well. This 
enables us to anticipate the behavior of the loss probability as n — * oo 
in some cases. 

Specifically, in the case g < 1 and n large, the loss probability is equal 
to conditional loss probability, provided that upon arrival at least m — 1 
servers are busy, multiplied by some constant. That is, as n large, the 
both abovementioned loss probabilities, conditional and unconditional, 
are of the same order. Following Abramov [2J (see also Choi et al [9]) 
this order is O(cr^). The precise result is given by Theorem 13. II below. 

If n large and g > 1, then the probability that arriving customer 
meets less than m — 1 customers in the system is small, and therefore, 
the loss probability should be approximately the same as the condi- 
tional stationary probability that upon arrival of a customer at least 
m — 1 servers are busy. That is, following Abramov [2] (see also Choi et 
al [9]) one can expect, that in the case of g > 1, the limiting stationary 
loss probability of the GI /M/m/n queue, as n ^ oo, should be equal 
to {g — l)/g for all m. Can the last property of asymptotic indepen- 
dence of m be extended as n increases to infinity and g approaches 1? 
The paper provides the condition for this asymptotic independence in 
this case. 
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3. Formulation of the main results 



Theorem 3.1. If g > 1 then for any m > 1 



lim pm,n = - — -■ (3.1) 

n^oo g 

If g = 1 and g2 = J^{fixYdA{x) < oo then for any m > 1 

lim npm,n = TT- (3-2) 

n— >oo Z 

If g = 1 and g^ = {fixY dA{x) < oo then for large n and any m > 1 
we have 

P^,n = ^ + 0[^). (3.3) 



If g < 1 then for Pm,n we have the limiting relation: 

lim ^^^^ = Kjn[i- + jjma'djm — j^mam)], (3.4) 

n—*oo (T" 

where «'(■) denotes the derivative of a{-), 



K 



^ 1 - gjj m l - a™ - 7J 



^ (1 - ifj) m(l - a„) - J 

/■oo 

= / e-'^^'^dAix), 







1 -^j 



anc? am is the least in absolute value root of functional equation: 

z = a{fim — fimz). 

Theorem 13.11 shows that if g > 1 then the limiting stationary loss 
probability is independent of parameter m. If g = 1 and g2 < 00 then 
lim„^oo npm,n is independent of parameter m as well. The proof of (13.11) 
seems to be given by simple straightforward arguments (extended ver- 
sion of the heuristic arguments of Section 2). Nevertheless, all results 
are proved by reduction to the abovementioned theorems of Takacs 
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[2lj and Postnikov [l7j given in the Appendix. The most significant 
result of Theorem 13.11 is (13.41) . This result is then used to prove the 
statements of Theorem 13.21 on the behavior of the loss probability as 
the load approaches 1 from the left. 

This behavior of the loss probability is given by the following theo- 
rem. 

Theorem 3.2. Let g = 1 — e, where e > 0, and en C as n oo 

and e 0. Assume that = g3{n) is a bounded sequence in n, and 
there exists g2 = lim^^oo Q2{n). In the case where C > for any m > 1 
we have 

Pm,n=Y^^Z^[l+0il)]. (3.6) 

In the case where C = for any m > 1 we have 

Pm,n = ^ + o(-). (3.7) 

Theorem 13.21 shows that as g approaches 1 from the left, the loss 
probability Pm,n becomes independent of parameter m when n large, 
and the asymptotic behavior of the loss probability is exactly the same 
as for the GI j M jl/n queue. 

4. Derivation of the recurrence equations for the loss 

probability 

For the sake of convenience, in this section we keep in mind that the 
first m — 1 states of the GI / M /m/n queue-length process form one 
special class. If an arriving customer occupies one of servers, then the 
system is assumed to be in this class, and the states of this class are 
numbered 1, 2,..., m. Otherwise, the system is in the other class with 
states m + 1, m + 2,..., m + n, where the last state, m + n, is associated 
with a loss of an arriving customer. 
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For example, if n = 0, then the second class of the GJ/M/m/0 
queueing system consists of one state only. 

Let us now build the recurrence relation similar to those of the 
GI/M/l/n queue. 

We start from the GJ/M/1/0 queue. For this queue we have 

7ri,o = — , (4.1) 

^0,1 

where 

ro,i = = / e-''MA(x). 
Equation (14. ip formally follows from the recurrence relation ttq^i = 
'"o,i7ri,o, where ttq,! = 1. The loss probability for the GJ/M/l/O queue 
is equal to 

1 

Pi.o = = Vi- 

7ri,o 

Before considering the case of the GI /M/m/n queue, notice that 
the value Tim-k,k has the meaning of the expected number of arrivals 
to the stationary system up to at the first time an arriving customer 
finds m — k servers busy and k remaining servers free. 

In the case of the GI/M/2/0 queue, by the total expectation formula 
we have 

TTl,! = ro,27r2,0 + n.lTTl,!, 

where 



oo 



= 2 / [1 -e-'^^]e-'^MA(a;), 



oo 



ro,2 = ^2= e-'^-dAix). 
Jo 

Then, by the total expectation formula, the recurrence relation for 
the GI/M/m/0 queue looks as follows: 

m—l 

TTm-l,! = ^^fc.m-fcTTm-fc.A:, (4.2) 

A;=0 
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where 



It is well-known that 




t(":)c. (4.3) 



and the loss probability is 



Pmfl — 



1 




(4.4) 



(see Cohen [TT], Palm [15], PoUaczek [16], Takacs [22] as well as Bharucha- 
Reid [5]). 

A relatively simple proof of (14 .Sp and (I4.4p can be found in Takacs 
|22j . It is based on another representation than (14.21) . For our further 
purposes, representation (14.21) is preferable. Representation (14.21) is a 
recurrence relation of the convolution type (II. 5p . and in the following 
it helps us to reduce the problem to the abovementioned combinatorial 
results of Takacs [23]. Once this is done, we apply then the Tauberian 
theorem of Postnikov [17] . 

Let us now consider the GI /M/m/n queueing system. In the case of 
this system with n > 1 we add an additional subscript to the notation. 
Specifically, rk,m-k,o = rk,m-k, and for j <n 



In addition, the value Tr^^on denotes the expected number of arrivals 
into the stationary GI/M/m/n queue up to the first loss. (Replacing 




and 
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the indexes n with j has the same meaning for the GI / M /m/ j queue.) 
Also, we will use the notation vr^^^ ^ g instead of the earlier notation 
T^m-k,k for the GI / M /m/Q queue. 

Then the recurrence relation associated with the nth series looks as 
follows: 

j m— 1 

1=0 k=l 

j = 0, l,...,n- 1, 

where 



(n) 

-i,i,0 

■ 



m—i,i,0 



rk,m-k^i+l,OTrm~k^i+l.k+i-l (4rm,0 - 1)' (4-6) 



1,2, ...,m - 1, 



and the second sum of fl4.5l) is equal to zero if m = 1. Moreover, if 
m = 1, then we do not longer need the upper index (n), showing the 
series number, and equation fl4.6l) . It is not difficult to see, that for the 
given series n, the recurrence relations (14. 5 p and (14.61) form a recurrence 
relation of the convolution type given by (II. 5p . In the next section we 
prove relations (13.ip - (]3.3p of Theorem 13. 1[ 

5. The proof of (ICT- flST^ 

First of all note, that 

m— 1 oo 

r 0,m,l 



1=0 k=l 1=0 

= V fe— ^^dA(a:) 
1=0 ^• 

= ^"X]e-"^^^^^:^^dA(x) = 1 
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Therefore, one can apply the theorem of Takacs pi] (see Appendix). 
Let 7i denote 

oo 

7i = ^lro,m,i- (5.1) 
1=1 

Then also 

n m—1 

lim \y lrQ^rn,i +y {n + k)rk ,m—k,n = 7i- (5-2) 

n— »oo L ' 

This is because 



1=1 k=l 



X] 



as n — >• OO. 

According to (15. ip and (15. 2p we have 71 = m/i/A, and therefore, 
7i = Then according to theorem of Takacs |23], given in the 

Appendix, in the case oi q > 1 we obtain 

n^oo i — 7l Q ~ i- 

Then, in this case of the limiting loss probability as n — > cxd we obtain 
lim pm,n = lim .\ = - — -. 

n^oo n— >oo ' Q 
m,n,0 

Similarly to ( 15. ip . Let 72 denote 

00 

72 = ^/(/-l)ro 

Then also 

n 

lim [ V /(/ - l)ro.m,i + y"{n + k){n + k-l) 



1=2 



m—1 



^k,m—k,n 

1=2 k=l 



72, 
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and < oo as 72 < oo. Indeed, as n — > oo, 

{n + k){n + k-l) 



X 



<{n + k){n + k-l) e-^irn-k)Am^^^^^^ ^ Q_ 

Therefore, in the case where g = 1 and Q2 = {fimx)'^dA{x) < 00 we 
obtain 

hm npm,n = TT- (5-3) 
Limiting relation (15. 3p can be improved with the aid of the Tauberian 
theorem of Postnikov [17] (see Appendix). In the case where g = 1 and 
§3 = {fimx)^dA{x) < 00, for large n we obtain 

Pm,n = — + ( — 



2n V n' 
Indeed, let 73 denote 

00 

73 = 5^^(^-l)(^-2)ro,m,z. 

1=3 

Then also 

n 

lim \\2l{l-l){l-2)ro,m,i 

n— >oo L ' 

m— 1 



+ ^(n + /c)(n + A; - l)(n + /c - 2)rfc 

,m—k,n 

k=l 

and ^3 < 00 as n — >■ 00. Indeed, as n — > 00, 

{n + k){n + k - l){n + k - 2) 
X (^"^^ j\-{m-k),.^ £ (^f^)!_!(e-M- _ e-'^-)^m/xdM}dA(i 

< {n + k){n + k ~ l){n + k ~ 2) \ 



m 
k 

"e-^(— ^•)-^^:^^dA(x)->0. 
n! 
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Thus, (E^D and (EJ} follow. 



6. The proof of (3.4) 

Whereas f l3.ip -( l373l) are proved by immediate reduction to the known 
results associated with (11.51) . the proof of (13. 4p requires special analysis. 
In order to simplify the analysis let us concentrate our attention to the 
constant in relation (13.41) . Multiplying this constant by (1 — am) 
we obtain 

(1 ^m) K^m 
-l-'^m ^ (1 - V^i) m(l - (Tm) - jJ 

The constant Km, given by (16. ip . is well-known from the theory of 
GI/M/m queueing system. Specifically, let pj, j = 0,1,..., be the 
stationary probabilities of the number of customers in this system im- 
mediately before arrival of a customer. It is known (e.g. Bharucha-Reid 
[5], Borovkov [6j) that for all j > m 

p, = Km<^. (6.2) 

Now, in order to prove (13. 4p let us write a new recurrence relation, 
alternative to (14.50 . For this purpose, join the first m states of the 
GI/M/m/n process to a single state and label it 0. Other states will be 
numbered 1, 2,..., n. In the new terms we have the following recurrence 
relations 

nf) = X:ro,^,nf_)^„ (6.3) 

with some initial value IIq"^ for the given series n. For example, for the 
series n = we have 



i=0 ^ ' 
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(see relation (14.31) ). A formal application of the theorem of Takacs [21] 
(see Appendix), applied to recurrence relation (16.31) . for large n yields: 

lim " = - . (6.4) 

n^oo j]in) 1 _|_ ^rna'{fim — ^mam) 

Let us now find lim^^oo nQ""*. Notice, that for j > m the probability 
Pj can be rewritten as follows. From (16. 2p we have: 

= K^a^~^^ = K^ai~-\1 - am) = (6.5) 

where Pj is the conditional probability for the GI/M/m queue, that an 
arriving customer finds j customers in the queue provided that upon 
arrival at least m — 1 servers are occupied. The conditional probability 
Pj-m coincides with the stationary queue-length distribution immedi- 
ately before arrival of a customer in the GI/M/1 queue given under 
the expected service time {fim)~^. Km is the stationary probability 
for the GI/M/m queue, that upon arrival at least m — 1 servers are 
occupied. 

From the theory of Markov chains associated with GI /M/l/n queue 
it is known (e.g. Choi and Kim [8j) that the j-state probability imme- 
diately before arrival of a customer is (7r„_j — 7r„_j_i)/7r„, where vr^ is 
given by (11.31) . and in turn the loss probability is determined by (II. ip . 

Then for the same j-state probability of GI/M/1 queue with the 
expected service time {fim)~^ we have 

tt('^) _ tt(") 

P, = lim (6.6) 

J-J-n 

In turn, by (16. 5p and (16. 6p . the j+m-state probability of the GI/M/ m 
queue is determined as follows: 

Pj+^-A„^iim 

J-J-n 
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lim n 



(n) 



(6.7) 



In view of (16.71) and (16.41) and according to Takacs theorem [2^ we 
obtain: 



lim 

n— »oo 



1 



qK„ 



Now, taking into consideration that the loss probability 

1 



Pm,r, 



J-J-n 



(6.8) 



we obtain statement (13.41) of the theorem. 



7. The proof of Theorem lOl 

It was shown in Subhankulov [21], p. 326, that if = 1 + e, e > 
and e 0, Qs^n) is a bounded sequence, and there exists §2 = 
lim„^oo Q2{n), then 

am = l-^ + 0{e'), (7.1) 

where = o'mln) is the minimum in absolute value root of the func- 
tional equation z = a{^m — fimz), \z\ < 1, and where the parameter 
/i and the function a{z), both or one of them, are assumed to depend 
on n. (Asymptotic representation (17. ip can be immediately obtained 
by expanding the equation z — a{fim — fimz) = for small z.) 
Then, after some algebra we have 

[1 + fima'{fim — jimam)] = £ + o{e), (7.2) 

and 

< = e-^^/^ni + o(l)]. (7.3) 
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In view of (17.11) . the term 



, 1 — cr„ 



^ (1 - (fj) m(l - (T„) - j 
has the order 0{e). Therefore, for term fl3.5p we have 

K^ = l + 0{e), (7.4) 

and in the case where C > 0, in view of (1 7. 2 1) -( 17^ and (16. 8p . we obtain: 

Pm,n = ^_^_,C/J ^ + 0il)]. (7.5) 

(13.61) is proved. 

The proof of (13.71) follows by expanding the main term of the asymp- 
totic expression of (17. 5p for small C. 
Theorem 13.21 is completely proved. 



8. Numerical example 

In this section a numerical example supporting the theory is provided. 
Specifically, we simulate D/M/l/n and D/M/2/n queues and check 
statements (13. 7p and (13. 6p of Theorem 13.21 numerically. The results of 
simulation are refiected in the table below. The value g is taken 0.999, 
so that e = 0.001 The value n varies from 10 to 50, and parameter 
C = en varies from 0.01 to 0.05 The theoretical values of the loss 
probability for these n are calculated by (13. 7p . There are also the loss 
probabilities for n = 100. The theoretical value for the loss probability 
related to this case is calculated by (13.60 . 

The table is structured as follows: Column 1 contains the values 
of parameter n, Column 2 contains the theoretical values for the loss 
probability given by (13. 7p and (13. 6p . Column 3 and 4 contain the loss 
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Loss probabihty 


Loss probability 


Loss probability 


n 


theoretical 


simulated for 


simulated for 






D/M/l/n queue 


D/M/2/n queue 


10 


0.0501 


0.0426 


0.0390 


15 


0.0334 


0.0292 


0.0275 


20 


0.0251 


0.0221 


0.0211 


25 


0.0200 


0.0180 


0.0173 


30 


0.0167 


0.0151 


0.0146 


35 


0.0143 


0.0128 


0.0124 


40 


0.0125 


0.0111 


0.0108 


45 


0.0111 


0.0098 


0.0096 


50 


0.0100 


0.0087 


0.0085 


100 


0.0045 


0.0040 


0.0039 


Table 1 . The comparison table of the 


OSS probabilities 



for D/M/l/n and D/M/2/n queues 



probabilities obtained by simulation for the D/M/l/n and D/M/2/n 
queueing systems respectively. 

As we can see from this table the difference between the loss prob- 
abilities of the single-server and two-server queueing systems obtained 
by simulation is not large, and difference between these loss probabil- 
ities decreases as n increases. As n increases the both simulated loss 
probabilities approach the theoretical loss probability. 
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APPENDIX 

In the appendix we recall the main results on asymptotic behavior 
of the sequence Q„, as n ^ oo (see relation (11.51) ). 

Denote f{z) = J2Zo 1^1 < 1> 7* = J2T=^ ( IlLj-^+i ^) fj- 

Theorem Al. (Takacs t24j, p. 22, 23.) //71 < 1 then 

lim Qn = , ^° . 

n— »oo 1 — 'Ji 

//71 = 1 and 72 < 00, then 

Qn 2(^0 

hm — = . 

n^oo n 72 

//71 > 1 then 

Qo 



lim 

n— >oo 



Qn 



5«(1 - f(5)) 



1-71 



where 6 is the least in absolute value root of equation z = f\z). 

Theorem A2. (Postnikov p!7], Section 25.) J/71 = 1 and 73 < 00, 

then as n ^ 00 

Q„ = ^n + 0(logn). 

72 
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